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ABSTRACT 

We study the growth of large scale structure in two recently proposed non-standard 
cosmological models: the brane induced gravity model of Dvali, Gabadadze and Por- 
rati (DGP) and the Cardassian models of Freese and Lewis. A general formalism for 
calculating the growth of fluctuations in models with a non-standard Friedman equa- 
tion and a normal continuity equation of energy density is discussed. Both linear and 
non-linear growth are studied, together with their observational signatures on higher 
order statistics and abundance of collapsed objects. In general, models which show 
similar cosmic acceleration at z ~ 1, can produce quite different normalization for 
large scale density fluctuations, i.e. ag, cluster abundance or higher order statistics, 
such as the normalized skewness S3, which is independent of the linear normalization. 
For example, for a flat universe with Q, m ~ 0.22, DGP and standard Cardassian cos- 
mologies predict about 2 and 3 times more clusters respectively than the standard A 
model at z = 1.5. When normalized to CMB fluctuations the as amplitude turns out 
to be lower by a few tens of a percent. We also find that, for a limited red-shift range, 
the linear growth rate can be faster in some models (eg modified polytropic Cardassian 
with q > 1) than in the Einstein-deSitter universe. The value of the skewness S3 is 
found to have up to ~ 10% percent variations (up or down) from model to model. 

Key words: cosmology: theory - cosmological parameters - large-scale structure of 
the Universe 



1 INTRODUCTION 

The suggestion that the universe is undergoing late time acceleration (eg Efstathiou, Sutherland & Maddox 1990; Riess et al. 
1998; Perlmutter et al. 1999) has provoked a number of cosmological models that can reproduce the appropriate evolution. 
The obvious solution is clearly to add the cosmological constant to the equations determining the evolution of the universe 
(see e.g. Weinberg 1989). This, as is very well known, is problematic from the point of view of particle physics due to the 
smallness of the required constant. It is then interesting to consider scenarios where instead of the cosmological constant, 
there is a new type of mechanism that can explain the acceleration of the universe. In particular, in this paper we consider 
two recent suggestions: the Dvali-Gabadadze-Poratti (DGP) brane induced gravity -model (Dvali, Gabadadze & Porrati 2000) 
and the Cardassian model(s) of Freese and Lewis (Freese & Lewis 2002). Both models lead to late-time acceleration with a 
non-standard Friedmann equation and with no explicit cosmological constant. The two scenarios are, however, fundamentally 
different in their properties, due to the fact that the DGP-scenario is a truly higher dimensional scenario while the Cardassian 
is an effective description. 

Both of these scenarios have been studied from the point of view of observational constraints (Deffayet et al. 2002; Wang 
2003; Gondolo & Freese 2002; Sen & Sen 2002, 2003; Zhu & Fujimoto 2002, 2003) mainly coming from the SNIa (Riess et 
al. 1998; Perlmutter et al. 1999) and CMB data. The constrains are not strong enough to exclude either of the models as 

* tuomas@ecm.ub.es 
f gazta@ieec.fcr.es 
I manera@ieec.fcr.es 



2 T. Multamaki, E. Gaztanaga and M. Manera 



an alternative to the standard A-cosmology. Are they compatible with other current observations, such as the large scale 
structure? 

In this paper we study the DGP and Cardassian scenarios from the point of view of large scale structure growth. Since 
these scenarios modify gravity on large scales, the growth is non-standard. The linear growth of perturbations in the original 
Cardassian scenario was also briefly considered in (Gondolo & Freese 2002). Our approach is much like that in (Gaztanaga 
& Lobo 2001), where the gravitational growth was studied in a number of non-standard scenarios. However, the Einstein's 
equations are modified on large scales, so we must alter our approach accordingly. 

This paper is organized as follows: In Section 2 we introduce the two scenarios and briefly consider the experimental 
constraints. In Section 3 we introduce the formalism of studying gravitational collapse in a matter dominated cosmology 
where energy density of matter is conserved but the Friedmann equation is arbitrary. The linear and non-linear growth of 
perturbations in a standard A-cosmology is also discussed in this Section. In Sections 4 and 5 we study gravitational collapse 
in the DGP and Cardassian scenarios and present the results. Observational constraints coming from the linear and non-linear 
aspects of structure formation are briefly considered in Section 6. The paper is concluded in Section 7. In the Appendix, the 
growth of perturbations in the original Cardassian model is studied analytically. 



2 THE TWO SCENARIOS 

2.1 Brane induced gravity model 

The brane induced gravity model (Dvali et al. 2000) offers an alternative explanation to the observed acceleration by a large 
scale modification of gravity due to the presence of an extra dimension (Deffayet et al. 2002a, b) . In this scenario there is hence 
no need for an explicit non-zero cosmological constant. At large enough scales gravity sees the full space-time, i.e. the brane 
and the bulk, and is therefore modified from the standard 1 /r 2 law 1 . 

In the DGP-scenario the Friedmann equation on the brane is (Deffayet et al. 2002a): 



where H = a/a, K = 0, ±1 is the curvature constant, a the scale factor of the universe and we have defined k 2 = (8ttG)/3, 
R = 2r c . The evolution of the scale factor is standard as long as the energy density of matter dominates i.e. k 2 p 3> l/(4r 2 ). 
Since p decreases with time due to the expansion of the universe, the r~ 2 -term becomes dominant at some point, and acting 
as an effective cosmological constant, it leads to an accelerating universe. At late times the scale factor grows exponentially, 
a ~ exp(t/r c ). The continuity equation is unchanged, 

p + 3H(p + p) =0. (2) 

We are interested in the growth of large scale structure which takes place in a matter dominated universe and hence we 
assume that p ^> p and therefore the continuity equation tells us that p ~ a~ 3 . We can view the DGP-scenario as standard 
cosmology with a perfect fluid with the normal properties. The only difference is the non-standard Friedmann equation. 

We define the cosmological quantities as usual 

^ M = -JJT ( 3 ) 

= -Jm ( 4 ) 



R 2 H 2 
K_ 
til 



A* = 772- (5) 

where -Ho is the current Hubble rate and po the current density of matter. Since we are interested in cosmology in the matter 
dominated universe, we can relate Qm to p by (we take ao = 1) 

_ K 2 p _ pa 3 



The Friedmann equation, Jy, can then be written as 



tfo 2 ((J^+^ + v^) 2 + ^). (7) 



1 It was later realized that in the DGP scenario gravity is modified also on small scales and may be detected in anomalous precession 
of orbiting bodies in the Solar System (Lue &; Starkman 2002; Dvali, Gruzinov & Zaldarriaga 2002) 
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Note that the normalization condition differs from the usual one, 

i = (yn^r+Ti^+vn^y + n K . (8) 

From now on we assume that we live in a flat universe, K = 0, so that Friedmann equation can be written as: 



The normalization condition simplifies to 

n M + = 1, (io) 

from where it is clear that in order to have Qm > 0, Qr must be restricted to the range 

0<fi fl <i (11) 

In (Deffayet et al. 2002a) this model was tested with SNIa and the CMB data. It was found that the model is in agreement 
with data 2 and the preferred parameter values for a flat universe are 

n M = o.is±g:gS. «h = o.i7tg;jg. (12) 

2.2 Cardassian models 

In the Cardassian 3 models (Freese & Lewis 2002; Gondolo & Freese 2002) the Friedmann equation has the general form 

H 2 =g(p M ), (13) 

where the pu is the energy density of ordinary matter and radiation. The universe is assumed to be flat and that there is 
no new type of matter nor a non-zero cosmological constant. The function g is assumed to approach the standard form, k 2 p 
at early times, including nucleosynthesis, and at late times, z < 0(1), to give accelerated expansion in accordance with the 
supernova observations ll9^ . Since the behavior of the function is different at different values of pM, there is an associated 
scale, pc, or red-shift, z eq , in the function g that determines when the evolution is standard and when the non-standard terms 
begin to dominate. The exact form of the model can hence vary, as long as it satisfies the aforementioned constraints. Some 
alternative forms of the Friedman equation are presented in (Freese 2002). 

As an example, consider the following (original Cardassian) form of g(p) (Freese & Lewis 2002) (we omit the subscript 
M from now on): 

H 2 = K 2 p + Bp n , n< |. (14) 

At early times the universe is dominated by the « 2 p-term (provided that B is small enough at the time of interest) and at 
late times the p n -term becomes significant, providing acceleration compared to the standard case. The term of the form p n in 
the Friedmann equation, and hence the Cardassian model(s), is motivated by considering our universe as a brane embedded 
in extra dimensions (Chung & Freese 2000) (this idea was, however, critically reviewed in (Cline & Vinet 2002)). In terms of 
the scale pc, Eq. 1141 can be written as 

H 2 = K 2 p[l + (— ) n_1 ], (15) 

PC 

hence B — K 2 pQ~ n . In a matter dominated universe this is conveniently parametrized by the red-shift at which the two terms 
are equal, z eq , 

H 2 = K 2 p[\ + (1 + *e,) 3Cl ~ B) (— ) n-1 ], (16) 

Po 

where po is the current energy density of matter. 

The requirement that n < | is actually more universal as one can see from considering the acceleration of the scale factor 
from the general Cardassian form, Eq. 11311 : 

I = ff(p) - \p9{p), (17) 

where it has been assumed that the energy density of ordinary matter is conserved and has no pressure i.e. p + 3Hp — 0. If we 
wish to have late time acceleration, a /a must be greater than zero at late times (when the assumption that P = becomes 
more and more exact) and hence at late times the inequality 



2 This result was, however, disputed by Avelino & Martins (2002) who argued that this scenario is already strongly disfavored. 

3 Humanoid-like race from Star Trek, see e.g. www.startrek.com 
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g(p)iate < Bp* (18) 

must hold. Therefore, in order to have late time acceleration, the Cardassian function, g(p), must grow more slowly than p3 
at late times. 

Recently, another Cardassian model has been studied in more detail with respect to the future SNIa observations (Wang 
et al. 2003). In the Modified Polytropic Cardassian (MPC) model, the Friedman equation is 

#Wp[i+(A) ]V ?i (19) 

where pc is again the energy density of matter at which the non-standard terms begin to dominate and q > is a parameter 4 . 
The original Cardassian model is a special case of the MPC model with q=l and hence we shall concentrate on this more 
general model in this paper. The Friedmann equation in the MPC-scenario in a matter dominated universe can equally well 
be written in terms of the red-shift at which the Cardassian terms start to dominate: 

H 2 = k 2 p(l + (1 + z eq f ^ « {f) in ~ 1)q ) ' = Hl^- (l + (1 + z eq f^"a^-^) ' , (20) 

where po is the energy density of matter today. 

In a flat universe, which we assume to the be case throughout this paper, the observed matter density of the universe, 
£Im b , can be related to z eq by (Wang et al. 2003) by 

i + Zeq = mt, a r q - i]V(s«(i-n))_ (21) 

The MPC-model is constrained by the supernova observations as well as the CMB (Wang et al. 2003). The experimentally 
allowed (n, q) parameter space is large and there is a degeneracy along the q axis at for q > 10 (when = 0.3) (Wang et 
al. 2003). 

The original Cardassian model has also been constrained in other works: in (Zhu & Fujimoto 2002) the angular size 
of compact radio sources at different red-shifts was considered and in (Sen & Sen 2002) by the SNIa, as well as in (Zhu & 
Fujimoto 2003), and CMB data. 



3 GRAVITATIONAL GROWTH 

Before studying the gravitational growth of structures in the two aforementioned scenarios, we first briefly introduce the usual 
formalism. 

Raychaudhuri's equation for a shear free fluid with four- velocity u M , is 
e + ie 2 = i? M „uV (22) 
where 

6 = V M n". (23) 
Choosing the coordinate system such that the four- velocity of the fluid, w M , is 

u M = (l,dx + v), (24) 
where v is the peculiar velocity. Hence, 

9 = 3- + -, (25) 
a a 

where we have defined 9 = V- v. Assuming that the geometry of the universe is of the standard FRW-form, the Raychaudhuri's 
equation simplifies to 

e + ^Q 2 = 3(H + H 2 ). (26) 

Note, that in order to study the growth of perturbations, on the LHS of Eq. 1261 . we are using the background quantities and 
on the RHS the local, perturbed quantities. Hence, using I25L we get 

- + -H+~=3(H + H 2 )-3(il + H 2 ). (27) 
a a 3 a z 



4 Incidentally, the Q-continuum is another race from Star Trek, existing in extra dimensions (www.startrek.com) 



5 



Up to this point, we have not assumed anything but that we live in a FRW universe filled with a perfect fluid and hence 
Eq. appl ies both to the DGP- and Cardassian-scenarios. In particular, we have not assumed any connection between the 
geometry and the energy content but simply rewritten the Raychaudhuri's equation in terms of the scale factor. This is useful 
since now we can study the evolution of the density perturbations armed with the Friedmann equation and the continuity 
equation, without having to concern us with the Einstein's equations. Obviously, these considerations make only sense on 
large scales where the effects of the non-standard cosmology can be seen. 

In order to have a deeper understanding of the different scenarios, we recall the growth of perturbations in the standard 
scenario with a non-zero cosmological constant. 



3.1 Standard scenario with A ^ 

In the standard scenario, we know the Einstein's equations so it is instructive to see how one arrives to the same result by 
using them and Eq. 12711 directly. 
From the Einstein's equations 

- ^g^TZ = -k 2 T m „ + A&,„ (28) 
and the energy-momentum tensor of an ideal fluid 

T M „ = Pg^ + (P + p)u„u v , (29) 
it can be verified that in the FRW-metric 

R^uV = -| K 2 (p + 3P) +A. (30) 
On the other hand, by using the Friedmann equation 

H 2 = K 2 p+± (31) 
and the continuity equation 

p + 3H{p + P) =0, (32) 
we see that 

3(H + H 2 ) = -^K 2 (p + 3P)+A (33) 

i.e. we arrive to the same equation using only the Friedmann and the continuity equations. The Raychaudhuri's equation in 
the standard case hence takes the form 



1 



- + -H+--=--^(p-p + 3(P-P)). (34) 
a a 3 or 2 

In a matter dominated universe, the continuity equation for a non-relativistic (p ^> p) fluid can be written as (Peebles 
1993) 

^- + (1 + 6)6 = 0, (35) 
ar 

where 

<5(r,x) = ^-l (36) 

is the local density contrast and we have switched to conformal time dt = adr. Using I35|l . I|36[l and l|34fl we see that in terms 
of the conformal time we get: 

d 2 S , -,d5 4 1 dS 2 3 2 , 2 

d^+^-siT^ =r {l+5)5pa > (37) 

where H. = d(ln a) /dr. Rescaling the time variable once more, rj = ln(a), we arrive at 

dT 2+(2 + if dn ~ 3 T+sW =2 K{1 + s) V' ( } 

Using the notations ®, and recalling that in the matter dominated regime it follows from the continuity equation, Eq. 
1321 1. that p ~ a -3 , the Friedmann equation can be written as 

H 2 = H 2 [^f + Q A ]. (39) 
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Using this in Eq. I38H gives us the well known result 



dr) 2 



2-1- 



Qi 



dS _ 4 _L_ (d5\ 2 = 3 , Q M 
dr) Sl + SKdr) 2 [ ' Q m +u 3 Qa' 



(40) 



2 Q.m + cl 3 Qa- 

In order to determine how a small perturbation grows with time at different orders in perturbation theory, we expand S 



«=i 



* = £ 



A fa) 



<>o. 



where So is the small perturbation (and the expansion parameter). Using the expansion, we get the linear equation 
3 fijvf . ; 3 Qm 



D{ + (2 - 



-Di = 0, 



2f2 M + a 3 n A 2Q M +a 3 Q A ' 
which in an Einstein-deSitter (EdS) universe (Qm = 1, = 0) has the well known solution 



D 1 (rj) = Cie 17 + c 2 e' 



-3»j/2 



(41) 



(42) 



(43) 



The solution to the linear equation in the general case with non-zero A, can be expressed in terms of the hyper-geometric 
function as 



Di 



ci- 



7,3/2 



+ C 2 2 Fl 



1 11 

3'T' 



The second order equation is 



D'i + (2-% 



Qi 



2 Q M + a 3 Q A 



Qa 3 



Da 



2 Q M + a 3 Q A 



3D 2 



= 0, 



(44) 



(45) 



where it is understood that the linear solution is substituted. Similarly one can recursively go on to arbitrary order. 

The second order equation determines how Gaussian initial conditions develop non-Gaussian features and can be related 
to t he skewness of the densi ty field at large scale. The (/-order moments of the fluctuating field are related to the perturbations 
by teernardeau et al. 20021) 

m q = (5"), (46) 
which in term can be related to the connected moments, or cumulants, £ q . The normalized skewness is given by l|Bernardeau et al. 20021 ) 



_ & m 3 

which can be written in terms of the first and second order perturbations. For example at leading order: 

m 3 = <<5 3 ) ~ (S't) + 3(S 2 5l) + ... 

For Gaussian perturbations (Sf) = 0, so that we get: 

'd! 



S3 = 3- 2 



(47) 



(48) 



(49) 



In an Einstein-deSitter universe this coefficient can be calculated exactly and it is S3 = 34/7 « 4.86. 

To illustrate the effect of the cosmological constant, we have solved the linear and second order equations numerically. 
The linear evolution and the evolution of S3 are shown in Figs and [5] 

We see that the effect on the linear growth factor is significant as perturbations grow less when the cosmological constant 
is large compared to the energy density of matter. The effect of the cosmological constant on linear growth is well represented 
by the other figure where future evolution is shown. With a non-zero Qa, Si freezes and structures stop growing in the future. 

The effect on S3 is very small, as expected, so that even when Qm = 0.3, AS3 < 1%. Such a small change is clearly out 
of reach of present day observations as is discussed in more detail in Section 6. 



3.2 Gravitational growth with a general Friedmann equation 

In a general case with the Raychaudhuri's equation given by Eq. 127L the equation governing large scale structure growth in 
the matter dominated universe can be derived following similar steps like in the standard case. The result is 



ti , (o , JL\M_±_l_ f d6_2_ l + 6/, T - T , tt2 , t -h;,-2 
73" 2 d n 3l + S [ dT> JJ 



2 (2 + =^--— -(-y = -3^[(H + H^-(H + H )), (50) 
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Figure 1. Linear growth with for different values of f2j\ 




where again r\ — ln(a). Obviously one cannot progress further unless the continuity equation is known so that H can be 
calculated. With the continuity equation one can then express p in terms of p and S using Eq. I36i . The resulting equation can 
be expressed in terms of cosmological quantities, SIm, Oa,..., by writing the Friedmann equation in terms of Qs and noting 
that from Eq. JHJ we get 

P=^ M . (51) 

We then have an equation determining the growth of density fluctuations, 8, in terms of Qs. We can expand the H + H 2 -term 
in terms of S and then the whole RHS of Eq. 1501 as 

3^-((H + H 2 )-(il + H 2 )\ =3(l+S)J2 c nS n - (52) 

H ' n=l 

Note that there is no constant term in the expansion. 

Expanding the perturbation according to Eq. 1 1 U the linear equation is then 

D'[ + (2 + z=o)D' 1 + 3ciDi = 0, (53) 
H 

the second order equation 
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D'i + (2 + - l{D' x f + 3 Cl Z> 2 + 6(ci + c 2 )D? = (54) 

H 3 



and further orders are easily found. For example, in the standard case we see that 

_ i K 2 p M i n M 

1 2 # 2 2n M + a 3 fi A V j 

q = 0, i = 2,3,.... (56) 

Using these, the linear and second order equations, Eqs 1421 and 1451 1 are easily reproduced. 

Note that these expressions make only sense on large scales where the non-standard evolution of the universe can be seen. 
On small scales the evolution obviously must be according to the standard Einstein's equations. 



4 DGP 

As we have seen, the two equations that determine how density perturbations grow are the Friedmann equation and the 
continuity equation. In the DGP- scenario these are: 

= p + 3H{p + p) (57) 



In the DGP-case we cannot use Eq. t\'M)t since the Einstein's equations on the brane are modified from the standard 
four- dimensional equations as is apparent from the non-standard Friedmann equation. Another way of seeing the same thing 
is to consider in the DGP scenario the quantity appearing on the RHS of the Raychaudhuri's equation, Eq. 1261 . 



H + H> = -^(1 + RV J p+ ^ )iP + P) +«P + 1 f 2 P + ^ + i- (59) 

This is different from the result that one gets in the standard (p + 3p), with the standard Einstein's equations. In 

the DGP-scenario the Raychaudhuri's equation hence needs to be calculated directly using the general approach described in 
the previous subsection. 

In the matter dominated DPG-scenario: 



±(H + H 2 )^ ^ ' , U.oi 



«m I V^r- \/H%-+n R ) +4a 3 n fl I + x /m + n R 



The coefficients c; are easily calculated and the first two are, again expressed in cosmological quantities: 



Cl 



C2 = 



2 (Q M + M R a 3 ) (Q M + 4 Q R a 3 ) + J Mr + &f a 3 (5 Mm + 8 SIr a 

4(f2 M + ttfla 3 ) 2 



16K/!¥+i/!!fl + ^ (fi A / + fijja 3 ) 3 



Finally, we need to calculate the term appearing in front of the <5'-term in the perturbation equation (1501) : 



Clearly, in the limit fl R = 0, all these expressions reduce to the corresponding expressions in the Einstein-deSitter case. 

We can now study the growth of perturbations numerically. The initial conditions are chosen such that at a = 10 -3 , the 
standard exponential solution, Di ~ exp(7;), is reached. In calculating the second order perturbation, initial conditions are 
chosen such that the standard solution, constant S3 = ^ w 4.86 is valid from the beginning. In Fig.|3the linear growth factor 
(and the linear growth factor normalized with the scale factor) for different values of Qr, with Qm then determined by Eq. 
IIOI I. are shown as a function of the scale factor. The non-linear growth is shown in Fig. 2] 

The general form of the linear growth factor is similar to the cosmological constant case but here the effect is even more 
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Figure 4. Non-linear growth for different values of Qr 



pronounced. For the preferred value Qr = 0.17 fovali et al. 2000) . the growth of linear fluctuations up to now is suppressed 
by a factor of 0.54 compared to the EdS-case. In the A-universe with Qa = 0.7, the suppression is 0.78. Hence, there will be 
significantly less structure growth on large scales in the DGP-scenario than in a A-cosmology. Furthermore, the suppression 
begins earlier in the DGP-scenario as is visible from Fig. 0] 

The second order perturbation, or S3, also grows differently from the A-cosmology. The value of S3 starts to grow at 
earlier times and varies more than in the standard A-universe. However, the variation is still observationally insignificant 
being at best less than one percent. 



5 CARDASSIAN MODELS 

In this section will consider the growth of gravitational instabilities in the Modified Polytropic Cardassian model described 
by Eq. 11911 . The relevant equations in a matter dominated universe are 

= p + 3Hp (62) 

H 2 = K 2 p(l + (l + Zaq f( 1 -^(^-^) lq (63) 
V p / 

From these all the relevant quantities can be calculated straightforwardly. 
The H + H 2 -term is 
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Figure 5. Linear growth for different values of n, q = 1 



1 1 + (1 + z^-^iSn - 2){0 n ~ 1)q 

2 l + (l + z e ,)3( 1 -")'?(^-)(n-i)g 



and the first two a coefficients 

1 1 + (4n - 2 + 3(n - l) 2 q)X + n(3n - 2)X 2 



Cl 



C2 = ^(l + X)- 3 (n-l) 



1 + X 



((l-n)«-l)(l + 3(n-l)g) 



+ \2 - An - (n - l)(9n - 5)g + 3(n - 1) 9 JX + n(2 - 3n)X 

where we have defined X = ((1 + z e5 )a) 3(1 - n),J in order to shorten the otherwise lengthy expressions. 
In the original Cardassian case with q=l, these expressions simplify to 



(64) 

(65) 

(66) 
(67) 



Cl 



(-2 



-| + n(l 



n)(l- 



\3(l-n) _3(l-n) 



H-Cl+^SP-^fflSC 1 -^ 

|n(n - 1)(1 - |n)(l + z e9 ) 3(1 - n) a 3(1 - n) 



1 + (1 + 2 eg ) a ( 1 -")a 3 ( 1 -") 

Again, as a check it is easy to see that with n = 0, we recover the standard coefficients 15511 . 

Finally, in order to study the growth of perturbations, we need the coefficient of the 5' term: 

Sq(l-n) 3.}(l-n) 



(68) 
(69) 



2 + 



H 

H 2 



3 1 +n(l + z e 



(70) 



2 1 + (1 + z eq f^ 1 - n )a^ 1 - n ) 

We can now calculate the growth of perturbations in the original Cardassian scenario. Initial conditions are chosen like 
in the DGP scenario. 

The solution to the linear equation in the original Cardassian model can be expressed in terms of the hypergeometric 
function. The growing part of the solution is found to be 



D 1 (x) = e x 2 F 1 



2 + 3n 11 - 6n 



3(1 — n)x 



n < 



(71) 



6(1 - n) ' 6(1 - n) ' 

For the general case, no such solution is found. 

We have plotted the linear and non-linear growth factors for n — 0,0.1, ...,0.6 and q = 1,2,5, 10 as a function of a in 
Figs 151101 The linear growth factor is presented both on a logarithmic scale and scaled by the scale factor. In plotting these 
figures, we have assumed that S7^J S = 0.3, which then sets the value of z eq according to Eq. (12 1 > . 

From the figures it is clear how the Cardassian scenario is fundamentally different from the DGP-model and a A-cosmology. 
Looking at the behavior of 81 with n > (since n = corresponds to a cosmological constant) the linear growth factor actually 
begins to decrease in the future i.e. areas of higher density begin to rarify. This is true regardless of the value of q. Near 
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Figure 6. Linear growth for different values of n, q = 2 
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Figure 7. Linear growth for different values of n, q = 5 
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Figure 8. Linear growth for different values of n, q = 10. 
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Figure 9. Non-linear growth for different values of n, <7 = 1, 2. 
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the turning point there is some interesting behavior as is visible from the plots of 8i/a. With larger values of q we see that 
before the growth factor starts to decrease, structures grow more rapidly compared to the standard case. This was not seen 
in neither of the other models studied in this paper. The value of Si /a at present is dependent on the value of the parameters 
but can again be as small as 0.5. 

In Fig. q = 1 and we are hence looking at the original Cardassian model. This is of interest since in this case we can 
understand some of the features by analytical considerations. For example, we see that the value of n that corresponds to the 
slowest growth is n fa 0.4. If one considers the equation determining linear growth in more detail, it can be seen that in the 
limit a — > oo, the slowest growth actually corresponds to n = | with Di(a) ~ a~3 as is shown in the Appendix. 

Looking at the figures depicting the change in S3, Figs I§1 and 1101 we see that in the Cardassian model the evolution of 
the value of S3 is strongly dependent on the value of q. The shape and the magnitude of variation changes with q so that with 
larger q one sees larger variations. 

In the original Cardassian scenario the values of S3 grow compared to the EdS scenario. The overall shape of the curves is 
similar to that of the standard A ^ universe described in Sec. 13. II However, here the scale of the change in S3 is much larger 
than in the standard A ^ universe. The growth of non-linearities seems to be fastest for n ss 0.3. Analytical considerations, 
see Appendix, show that in the large a limit, the value of n corresponding to fastest growth at large a is again n = |. The 
magnitude of change in S3 is approximately two percent with n 
A-cosmology, is out of reach of present day experiments. 

However, as q is increased, also the variation of S3 grows. For example, with q 



0.3, which, although larger than in the DGP-scenario or 
10, we see that S3 can vary more than 
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Figure 11. Value of <rg in non-standard cosmologies, with Qm = 0.22, relative to the value for the Q\ = 0.78, C2m = 0.22 standard 
cosmology with erg = 1. 

10 percent from its EdS value. Also with q > 2, the values of S3 undergo a change from values smaller than S^ dS to values 
larger than Sf dS . 

6 OBSERVATIONAL CONSTRAINTS 

As mentioned in the introduction, the above non-standard models have been tested mostly against observations of SNIa and 
CMB. The analysis presented in this paper allows for a further comparison with constrains coming from linear and non-linear 
aspects of structure formation, i.e: <rg normalization, cluster number density and skewness of the density field. The shape of the 
initial spectrum of fluctuations P(k) as a function of wavenumber k, could also depend strongly on cosmological parameters 
(see eg Efstathiou, Sutherland & Maddox 1990). But this requieres additional ingredients to our models: ie knowledge of 
the primordial spectrum and the matter content. Here we concentrate only in aspects that relate to the matter dominated 
evolution, assuming some given initial shape P(k). We briefly sketch how this comparison can be done and leave a more direct 
parameter estimation for future work. 

6.1 Amplitude normalization: as 

The amplitude of density fluctuations is commonly characterized by the linear value of the rms fluctuations on a sphere of 
8Mpc/h. Observationally its value seems to be of order unity at z = 0. Form the analysis in the previous sections it is straight 
forward to predict how different as should be for a given cosmological model. One way to do this is to fix the normalization of 
all models to be equal at high red-shifts and see how different they are at z = 0.5. Figure [Til illustrates this point. We show 
the value of as normalized to the A model as a function of red-shift for DGP and MPC (n = 0.6 and 3 = 1) cosmologies, all 
with same value of Q m = 0.22. 

Values of as for other parameters can readily be obtained by comparing Figure with Figure [3] and Figures 15181 Note 
how for MPC with q > 1 we can also get as > 1. 

Current observational constraints of as usually come with several different assumptions. In most cases a shape for the 

5 In principle, this is similar to a CMB normalization, but note that in general a CMB normalization also requires some additional 
information on the shape of the spectrum, which might vary from model to model. Thus, here we choose to illustrate how erg changes 
due to the linear growth from a given fixed shape of the spectra. 
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Figure 12. Bottom: Comoving volume element as a function of red-shift normalized to the EdS model. The short and long dashed line 
shows DGP and MPC predictions for Qm = 0.22 (On = 0.15) and h = 0.71. The continuous line corresponds to the A cosmology with 
Qm = 0.22 (Qa = 0.78) and h = 0.71. Top: Number density of predicted clusters per square degree with mass M > 2 X 10 14 M so ; ar . for 
the same models. All models are normalized to erg = 1 at z = 0. 



linear spectrum is assumed to get an estimation, for example from cluster abundances or normalization of CMB fluctuations. 
If the primordial spectrum of fl uctuations is not s cale invariance or a simple power-law, such extrapolations could yield 
misleading values of og (e.g. see tearriga et al. 200ll) and references therein). Large scale structure in local galaxy catalogues 
also yield values og ~ 1, but one needs to quantify the bias b, or how well the selected galaxies trace the underlying mass 
distribution (see e.g. Gaztanaga 1995) . Recent CMB results from the WMAP (Spergel et al. 2003) mission fit to the (power- 
law) ACDM cosmology find: ag = 0.9 ± 0.1 (ag = 0.84 ± 0.04 with a running index, iSpergel et al. 20031) ). This value has 
been extrapolated from z ~ 10 3 to z — by using the linear growth fac tor Di(z) in Eq. 14411 . Significantly smaller values 
(ag ~ 0.6 — 0.7) have been found at low red-shifts by weak lensing (see e.g. l)jarvis et al. 200*21) ) and velocity fields (e.g. Willick 
& Strauss 1998). This apparent discrepancy, if real, could be easily accounted for by using the DGP or MPC cosmological 
models. 



6.2 The number counts of galaxies 

Press & Schechter formalism (Press & Schechter 1974) and its extensions (see e.g. Bond et al. 1991; Lacey & Cole 1993) predict 
the evolution of the mass function of collapse objects. These predictions are based on assuming Gaussian initial conditions 
and the spherical collapse model (which is closely related to the shear-free approximation used in Eq. I|22fl ). Despite the 
apparent limitations of these assumptions, comparison with realistic simulations show good agreement (see a recent review 
by Cooray & Sheth (2002)). In the standard Press-Schechter formalism, the comoving number density of collapsed objects 
(halos or clusters) of mass M is 

,»^,»,r l~2 f8c\ diner ( 8 2 c \pdM 

n(M)dM = - A /- — - — exp c - \'- r -^- (72) 

v 1 \ ty \ a J din M \ 2a 2 J M 2 x ' 

where a = a(R,z) is the linear rms fluctuation at the scale R corresponding to the mass M = i/3ivR 3 p, and p is the mean 
background. The value of S c corresponds to the value of the linear over-density at the time of collapse, i.e. when the non-linear 
over-density becomes very large S — > oo. This value can be found by solving Ea. l|37|l . For the standard Einstein-de Sitter case 
we have S c ~ 1.686. We find little difference in S c for the non-standard cosmologies, so for simplicity we will use the EdS value 
in all cases. 
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To translate the above predictions into observable quantities, such as the number density of clusters per unit red-shift, 
above a certain mass, we need to integrate over the comoving volume element dV/dz, which is also a strong function of the 
cosmology. Bottom panel of Figure \TI\ shows the comoving volume element (normalized to EdS case) in the DGP (dashed 
line), MPC (long dashed line) and A cosmology (continuous line). As can be seen in the figure, the A cosmology has about 
4 times more comoving volume by z ~ 1.5 than the EdS case, while DGP and MPC are only 3 and 2 times larger. Despite 
the smaller volumes, DGP and MPC predicts 2 and 4 times more clusters at z ~ 1.5 than the A cosmology because of the 
stronger freeze in the linear growth factor, which can be seem by comparing Figure^to Figure^! 

Note that there are four important factors in the PS predictions: i) Qm, which relates a given cluster mass M to R in 
a(R) value, ii) the shape of the power spectrum, which determines the a(R) curve, iii) the volume element, and iv) the growth 
factor Di(z) which gives the red-shift variation of a. In our analysis we fixed i) and ii) to the A cosmology with Qm — 0.22. 
Thus all differences in Figure WI\ are due to differences in volume and growth factors, which mark the distinction between 
standard and non-standard cosmologies. 



6.3 The skewness S3 

We have already shown the predictions for the normalized skewness, S3 in Eq. 1471 . in Figures |2"1 HI l9land !10l ("see also Eq. <A8t 
in the appendix). These are the unsmoothed values of the skewness at a single point. In practice, to relate to observations, 
one needs to take into account smoothing effects. For a top-hat window of a given radius R, smoothing results in a simple 
linear correction that is given by the slope 7 of th e variance smoothed at a radius R (Juskiewicz, Bouchet & Colombi 1993). 
In the standard cosmology: S3 = 34/7 + 7 (eg see l)Bernardeau et al. 2002h and references therein). The smoothing correction 
should be the same for non-standard cosmologies, as long as we work within the spherical collapse model (ie the shear-free 
approximation in Eq. 122H ). which gives the exact leading order contribution to S3 for a top hat window (see e.g. Fosalba 
& Gaztanaga 1998; Bernardeau et al. 2002). Thus, the difference between the standard predictions and the predictions from 
non-standard cosmologies should not be affected by these smoothing effects. Besides smoothing effects, S3 is also affected by 
systematic uncertainties of biasing (see §7. 1 in Bernardeau et al. 20 02). 

Current estimations for S3 (see §8 in ijBernardeau et al. 20021 ) ) agree with the standard predictions but with large un- 
certainties, of order 20% — 30%. Thus current observations on S3 can not be used to separate the different models. Next 
generation catalogues, such as the Sloan Digital Sky Survey, are expected to reduce the sampling variance uncertainties to 
5% (see e.g. Szapudi, Colombi & Bernardeau 1999; Verde et al. 2000). Thus, after removing biasing uncertainties, future 
measurements on S3 could be used to constraint some of the MPC parameters (e.g. see Figs. l9l and 1101 . 



7 CONCLUSIONS 

In this paper we have studied the growth of large scale structure in non-standard cosmological scenarios that can undergo 
acceleration without having an explicit cosmological constant. In order to do this, we have extended the standard spherical 
collapse formalism to account for the non-standard equations. All that is needed is the Friedmann equation and the assumption 
that the energy density of matter is conserved. The applicability of our approach is limited to very large scales, where Einstein's 
equations can be modified. Of the two studied scenarios, the DGP-model is more attractive since it is based on physics coming 
from brane cosmology whereas the MPC-model is a more phenomenological approach. In a flat universe the DGP-model has 
only a single free parameter, Qm, but this model is still found to be in accordance with the current observations on large 
scale structure. The MPC model has two additional parameters so that it is less surprising that observations on large scale 
structure do not rule out such a model. 

In the DGP- and MPC- models it was found that linear growth is inhibited due to the accelerated expansion, like in the 
standard A-cosmology. Depending on the values of the parameters, the suppression at present compared to the EdS-case can 
be as much as 50%. The non- linear growth was studied by considering the value of the skewness. For DGP, the deviation from 
the EdS-value is larger than in the A-case but still only of the order of one percent. In the MPC-scenario, however, the story 
is quite different and the value of skewness can vary up to 10 percent. 

In order to evaluate the significance of the non-standard evolution of large scale perturbations, one obviously needs to 
relate the predictions of the linear and non-linear growth rate to observations. As was discussed in Section 6, the linear 
normalization of large scale density fields can vary greatly. The value of erg can be much smaller in the MPC- and DGP- 
scenarios when compared to the standard A-cosmology with the same normalization at high redshift. Intriguingly, there is 
already some observational indication of a small value of as at low red-shifts but more observational data is needed in order to 
determine whether this effect is real. We also find that in the MPC-scenario, the linear growth can actually be faster than in 
the EdS universe for a limited range of red-shift. This is an interesting property since models that have late time acceleration 
typically lead to less linear growth at all times. In any case, observations of as at low red-shifts, combined with the precision 
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CMB data coming from the WMAP (and in the future Planck) mission, will be a powerful probe of non-standard cosmological 
evolution. 

Another quantity that is strongly affected by the non-standard linear evolution is the number counts of objects. As it 
was shown in Section 6, using the Press-Schechter formalism gives strong predictions for the number of clusters that differ 
significantly from one expects in a A-cosmology. Again, studying the number counts of clusters at different red-shifts gives 
another way to constrain non-standard cosmological models such as the DGP- and MPC-scenarios. 

The higher order statistics can also be modified by the non-standard effects. In this paper we have considered the 
normalized skewness, S3, whose evolution can be straightforwardly extracted in each model. The skewness gives another 
probe of possible atypical cosm ological evolution. Typically, th e effect on the skewness of different cosmological scenarios, e.g. 
A-cosmology and quintessence feenabed fc Bernardeau 200 lh . is very small and totally unobservable. However, as we have 
seen here (see also Gaztanaga & Lobo 2001), the effect can be much larger and can possibly be detected in the near future 
experiments such as the Sloan Digital Sky Survey. 

The growth of large scale structure is tightly bound to the overall evolution of the universe and hence to the particular 
cosmological model. As such, along with the SNIa and CMB data, it is a powerful tool to probe the space of possible 
cosmologies. Current large scale structure data agrees well with the two proposals that explain late time acceleration without 
a cosmological constant. This agreement is not trivial in the case of DGP, as there are no additional free parameters, once we 
fit to the SNIa data or the value of Q.m- We have shown here how upcoming surveys of large scale structure can be used to 
further constrain these non-standard models and differentiate then from the now standard A model. 
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APPENDIX A: GROWTH OF LINEAR AND SECOND ORDER PERTURBATIONS IN THE 
CARDASSIAN MODEL 

In this appendix, we consider the growth of linear and second order fluctuations in the original Cardassian model, q = 1, by 
analytical means. 

Al Linear growth 

The equation determining linear growth, Eq. 1531 . in the Cardassian model in a scaled form is 
<fD,( S l + ne^-J' Wfli . 3 2n(l-|n)e 3 ( 1 -^-l 



da; 2 V 2 1 + e 3 ^-")* J dx 2 e 3(i-n)s + 1 
where a; = ln(a/ac) = V — rjc- 

It is now interesting to look at the large \x\ limits. Since we are interested in the range of values where < n < |, it is 
clear that at, depending on the sign of x, the exponential terms in Eq. ijAlfl will either dominate or be negligible in the large 
x limit. When x is large and negative, i.e. when a/ac >C 1, Eq. ijAl^ takes the form 

i.e. the standard form, which has the usual solution 

£>r0) = A 1 e x + A 2 e-§ x . (A3) 

In the large positive x-limit, i.e. when a/ac ^> 1, it is obvious that the exponential terms will dominate and so that Eq. 
(| A ip> can be written as 

d^Di /n 3 . dD\ n 3 * _ „ / * .\ 

^ + (2 ~2 rl) ^ +3n(1 -2 n)Z?1 = a (A4) 
The solution to this equation is easily found and reads as 

Df(x) = B ie -? nx + B 2 e (3n - 2)x , (A5) 

where Bi are constants. Looking at the solutions, we see that with n = 1, the two solutions agree as expected. If n < |, the 
linear growth of perturbations will at some point stop and perturbations start to shrink. 

The slowest growth rate at large a/ac is easily deduced from Df(x) and occurs when — |n = 3n — 2 i.e. when n = | and 
hence D~[ (x) ~ e z x . 

A2 Second order perturbations 

The equation determining the growth of second order perturbations in the large x limit, in the negative large x limit we 
obviously reproduce the standard scenario again, is from Eq. I|54|l . 

d 2 L> 2 3 ,dD 2 3 8 fdD 1 2 3 2 

+ ( 2_ o"-)-5 \-3n{l - -n)D 2 - - 1 -j— ) + 3n(n + 1)(1 - -n)^ = 0. (A6) 



dx 2 v 2 ' da- ' 2 J 3\ dx ' y " 2 

The part that will dominate the linear solution Di depends on n, if n < |, Di ~ exp(— |na;), where as if | < n < 1, 
D\ ~ exp((3n — 2)a;). Let us first assume that n < |, in which case by substituting the appropriate solution of D\ into Eq. 
\A6l and solving the resulting equation, we get 

D 2 {x) = i(2 + n) e - 3na; + Cie^i"- + C 2 e {3n ~ 2)x . (A7) 
Hence, since n < |, we see that in the large a; limit, D 2 ~ e - ^™ 1 . Therefore, we expect that in this limit, 

S 3 ~ e* M . (A8) 

If we instead look at the values of | < n < 1, we see that S3 ~ e' 2-3 "' 21 , and hence the value of n that gives the largest effect 
rate of change of S3 is, again, at n = | with S3 ~ e^^. 



